precision='double';
for n=2:2:16
    A=rand(n,n);
    
    tic;
    [pp1a pp2a]=F(A);
    Ba=adjF(pp1a,pp2a);
    t1 = toc;

    tic;
    Bb=adjFFUsing6Loops(A);
    t2=toc;

    tic;
    Bc=adjFFUsing4LoopsAndH(A);
    t3=toc;

    reportTestResult(n,(norm(Ba(:)-Bb(:))./norm(Ba(:))),t1/t2,precision);
    reportTestResult(n,(norm(Ba(:)-Bc(:))./norm(Ba(:))),t1/t3,precision); 
end